Size, oxidation, and strain in small Si/Si02 nanocrystals. 
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The structural, electronic and optical properties of Si nanocrystals of different size and shape, 
passivated with hydrogens, OH groups, or embedded in a Si02 matrix are studied. The comparison 
between the embedded and free, suspended nanocrystals shows that the silica matrix produces a 
strain on the embedded NCs, that contributes to determine the band gap value. By including the 
strain on the hydroxided nanocrystals we are able to reproduce the electronic and optical properties 
of the full Si/SiC>2 systems. Moreover we found that, while the quantum confinement dominates in 
the hydrogenated nanocrystals of all sizes, the behaviour of hydroxided and embedded nanocrystals 
strongly depends on the interface oxidation degree, in particular for diameters below 2 nm. Here, the 
proportion of NC atoms at the Si/SiC>2 interface becomes relevant, producing surface-related states 
that may affect the quantum confinement appearing as inner band gap states and then drastically 
changing the optical response of the system. 

PACS numbers: 73.22-f; 71.24.+q,73.20.at; 78.67.Bf. 
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I. INTRODUCTION 

Silicon nanocrystals (Si-NCs) have known very in- 
tensive research activity since last decade because of 
their very promising light emitting properties. 1,2 The 
confinement of carriers in silicon structures leads to 
a strong enhancement of the radiative transition yield 
and to increase of the emitted photon energy. Em- 
bedding Si-NCs in wide band gap insulators is one 
way to obtain a strong quantum confinement (QC). Si- 
NCs embedded in a silica matrix have been obtained 
with several techniques as ion implantation, 3 ' 4 chemical 
vapour deposition, 5-8 laser pyrolysis, 9,10 electron beam 
lithography, 11 sputtering, 12,13 and some others. Alter- 
native preparation methods produces NCs with different 
surface terminations (D, F, C, and others) that earned 
some interest in the last years. 14-18 Moreover, also the 
role of nitrogen has been recently explored 19 . 

If the response of single NCs can be described as 
atomic-like, 9,11,20 in experimental samples no two NCs 
are the same, and collective properties have to be mea- 
sured. Experimentally, several factors contribute to make 
the interpretation of measurements on these systems a 
difficult task. For instance, samples show a strong dis- 
persion in the NC size, that is difficult to be determined. 
In this case it is possible that the observed quantity does 
not correspond exactly to the mean size but instead to 
the most responsive NCs. 21 Again, NCs synthesized by 
using different techniques often show different properties 
in size, shape and in the interface structure. The charac- 
terization of the effects of these quantities (size, interface 
configuration and passivation types) on the final proper- 
ties of the system, is of fundamental importance in order 
to determine the best condition for light emission and 
optical gain. 



Previous works 22-24 already observed that, while the 
electronic gap presents an oscillating behaviour with the 
size of the NC, the absorption threshold follows the QC 
for all the size considered. The effect is that the expected 
photoluminescence (PL) energy (that is associated to the 
electronic band gap) and the observed Stokes-shift (that 
is related to the difference between the absorption and 
emission energies) , cannot be simply described as a func- 
tion of the bare NC size, but the specific configuration of 
the interface should be taken into account. 10 
It has also been demonstrated that, for these systems, 
while the inclusion of local fields provide an important 
contribution to the absorption spectra, the excitonic 
and Coulomb corrections almost exactly cancel out each 
other. 22 This allows a simplification of the algorithms 
adopted in the calculations, still producing reliable re- 
sults. 

In this work we present first-principles calculations of the 
structural, electronic, and optical properties of Si-NCs of 
different size and shape, passivated with hydrogens, OH 
groups, and embedded in a Si02 matrix. The effects 
generated by the surrounding matrix on the electronic 
properties of the NCs (of different size) are analyzed. 
In particular, we compare a set of suspended (free) and 
embedded NCs in order to investigate the effects of the 
strain and the role of oxidation on the final properties of 
the system. 

The paper is organized as follows. A description of the 
theoretical methods and of the systems under investiga- 
tion is given in section II. The structural, electronic and 
optical properties for the embedded and suspended NCs 
are analyzed in section III. In section IV the effects of 
the NCs size and oxidation are discussed. Conclusions 
are presented in section V. 



FIG. 1: (color online) The final optimized structure of the Si32 in /3-cristobalite matrix (left panel), Si32-OH (center panel), 
and Si;!2-H (right panel). Red (dark gray) sticks represent the O atoms, cyan (gray) sticks represent the Si of the matrix, and 
the yellow (gray) thick sticks represent the Si atoms of the NC. 



II. THE METHOD 

The /3-cristobalite (BC) SiC>2 is well known to give 
rise to one of the simplest Si/Si02 interface because of 
its diamond-like structure. 25 The crystalline embedded 
structures have been obtained from a BC cubic matrix 
by removing all the oxygens included in a cutoff-sphere, 
whose radius determines the size of the NC. By center- 
ing the cutoff-sphere on one silicon or in an interstitial 
position it is possible to obtain structures with different 
symmetries. To guarantee a proper shielding of the in- 
troduced strain by the surrounding silica we preserved a 
separation of about 1 nm between the NCs replica. 26-28 
The result is a structure with the silicon atoms bonded 
together to form a crystalline skeleton with T^ local sym- 
metry before relaxation. In such core, Si atoms show a 
larger bond length (3.1 A) with respect to that of the Si 
bulk structure (2.35 A). No defects (dangling bonds) are 
present, and all the O atoms at the Si/Si02 interface are 
single bonded with the Si atoms of the NC. 
The optimized structure has been achieved by relaxing 
the total volume of the cell (see example in Fig. 1, left 
panel). The relaxation of all the structures have been 
performed using the SIESTA code 29,30 and Troullier- 
Martins pseudopotentials with non-linear core correc- 
tions. A cutoff of 150 Ry on the density and no addi- 
tional external pressure or stress were applied. Atomic 
positions and cell parameters have been left totally free 
to move. 

The hydroxided-strained NCs (s-Si^^-OH) are obtained 
by extracting the NCs together with the first interface 
oxygens from the relaxed NC-silica complexes, and then 
passivating the surface with hydrogen atoms (Fig. 1, cen- 
ter panel) . Also the case of hydrogenated-strained NCs is 
considered (s-SL^-H) , by replacing the OH groups with 
hydrogens (Fig. 1, right panel). In the last two cases, 
in order to preserve the strain, only the hydrogen atoms 



have been relaxed. The equivalent hydroxidcd-rclaxed (r- 
Siaa-OH) and hydrogenated-relaxed (r-SL^-H) NCs are 
obtained by a further optimization and relaxation of the 
whole systems. 

The goal is to distinguish between the properties that 
depend only on the NC from those that are instead in- 
fluenced by the presence of the matrix. The comparison 
of the results relative to different passivation regimes (H 
or OH groups) could give some insight on the role played 
by the interface region. Electronic and optical proper- 
ties of the relaxed structures have been obtained in the 
framework of the density functional theory (DFT), using 
the ESPRESSO package. 31 Calculations have been per- 
formed using norm-conserving pseudopotentials within 
the LDA approximation with a Ceperley- Alder exchange- 
correlation potential, as parametrized by Perdew-Zungcr. 
An energy cutoff of 60 Ry on the plane wave basis have 
been considered. 



III. EMBEDDED VS SUSPENDED NCS 

In this section we analyze the role of the surrounding 
silica matrix on the electronic and optical properties of 
the Si/Si02 systems. We compared the results of the 
composite Si/Si02 systems and the suspended ones, in- 
cluding the effects of the strain induced by the silica ma- 
trix. 

Four embedded system have been produced: the Siio 
and Sii7 structures have been obtained from a BC-2x2x2 
supercell (192 atoms), while for the Si32 and Si35 NCs 
the larger BC-3x3x3 supercell (648 atoms) has been 
used. 

Being the bonds near the interface much more strained 
with respect to those in the NC-core, 28 ' 32,33 the average 
strain tends to zero for large clusters. Therefore we 
adopted the most strained bond value (MSB) (calculated 



2 




2 4 6 8 10 12 

Energy (eV) 



FIG. 2: (color online) DFT-RPA imaginary part of the dielec- 
tric function for the Si32/Si02 (solid), S-S132-OH (dashed), 
r-Si 32 -OH (dotted). 



from the difference of the largest Si-Si distance of the 
NC, with respect to the bulk value), as the parameter 
representing the NC strain. For the four structures 
enumerated above we obtained MSB values ranging 
between 0.1 and 0.5 A. The corresponding suspended, 
hydroxidcd and hydrogcnated, strained and relaxed NCs 
were obtained following the procedure described in the 
previous section. 

The calculated energy differences between the highcst- 
occupied-molecular-orbital (HOMO) and the lowest- 
unoccupied-molecular-orbital (LUMO) for all the 
systems are presented in Table I. First of all, we note 
that the gap of the embedded NCs presents an oscillat- 
ing behaviour with the NC size (disregarding the QC 
rule), as already evidenced in other works. 22,34 For the 
hydroxidcd systems we note that, while the removal of 
the embedding medium produces a constant reduction of 
the gap of about 0.15-0.19 eV, the removal of the strain 
results in a gap enlargement that strongly depends on 
the system considered, ranging from 0.07 eV (Si35 NC) 
to 1.29 eV (Sii 7 NC). This result reveals that the NCs 
can be subjected to amounts of strain that differs very 
much one from another. Besides, for strained NCs, also 
the substitution of -OH terminations with -H produces 
gap enlargements ranging from 0.3 eV (Sin NC) to 3.06 





Siio 


Sil7 


S132 


Si35 


embedded 


1.77 


2.36 


2.62 


1.62 


hydroxided-strained 


1.60 


2.21 


2.45 


1.43 


hydroxided-relaxed 


1.97 


3.50 


2.93 


1.50 


hydrogenated-strained 


4.66 


2.51 


2.90 


3.14 


hydrogenated-relaxed 


4.61 


4.06 


3.54 


3.43 



TABLE I: HOMO-LUMO gaps for the embedded and sus- 
pended set of NCs. Units are in eV. 



eV (Siio NC). These strong variations of the gap due to 
strain removal and termination type will become clearer 
in the following, when the relationship between strain 
and oxidation will be taken into account. 
As already observed, 22 even for the completely relaxed 
hydroxidcd systems we have a strong discontinuity of 
the gap values with the NC size. Instead, the removal 
of the strain in the hydrogcnated systems produces a 
trend of the HOMO-LUMO gaps that follows the QC, 
as expected. 35 

To analyze in detail the role of the embedding matrix, 
as for the contribution to the absorption spectrum, like 
as for the induced strain, we consider the reference case 
of the Si 32 NC. The absorption spectra [expressed by 
the imaginary part of the dielectric function, calculated 
within the DFT-random-phasc approximation (RPA)] of 
Si32/Si02, S-S132-OH, and optimized r-Si32-OH are pre- 
sented in Figure 2. We note that the hydroxided-strained 
NC is able to reproduce very well the spectrum of the full 
Si/Si02 system in the 0-7 eV region, that is indeed as- 
sociated to the NC+interface contribution. 22,34 Instead, 
the removal of the strain produces an enlargement of the 
HOMO-LUMO gap, 36,37 and a consequent blue-shift of 
the absorption spectra in this region. For higher ener- 
gies (above 7 eV), the spectra of the suspended passi- 
vated clusters do not match that of the composite sys- 
tem, due to the absence of the silica matrix. These are 
very general results, that we have verified in all the cases 
of Si/Si0 2 NCs. As already demonstrated, 22,35,38 - 40 the 
substitution of hydroxyl groups with hydrogens results in 
a strong blue-shift of the absorption threshold, and the 
disappearance of the low-energy features that are related 
to the Si/Si02 interface. 

To conclude this section, we observe that the considered 
NCs present different configurations of the Si/Si02 inter- 



G 




0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.55 
MSB (A) 



FIG. 3: (color online) Dependence on the MSB value 
for (squares) HOMO-LUMO gap difference between the 
hydroxided-relaxed and -strained NCs, and (circles) oxidation 
degree. The lines are drawn to guide the eye. 
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FIG. 4: (color online) Stick and ball pictures of the final optimized set of hydroxided structures. Red (dark gray), white (light 
gray), and yellow (gray) spheres represent respectively the O, H, and Si atoms. 



face. In particular, the silicon atoms at the interface have 
a different number of nearby silica oxygens. Intuitively, 
we expect that a greater number of oxygen atoms could 
exert a larger strain on the interface atoms. By defining 
S7 as the number of average oxygens per silicon at the 
interface, i.e. as the ratio between the number of oxygen 
atoms and the number of silicon atoms bonded to them, 
we could compare the amount of strain with the oxidation 
fL By looking at Figure 3 we observe, by comparing the 
properties of the hydroxided-relaxed and -strained NCs, 
that both the oxidation degree (upper curve) and the 
HOMO-LUMO gap shift (lower curve) generally increase 
with the strain, as expected. A discrepancy is present for 
the Siio NC, that could be reasonably addressed to the 
fact that this is a borderline case in which all the Si atoms 
of the NC are localized at the interface. Also, the discon- 
tinuity could arise from the contribution of the angular 
strain (not considered here), and from the non uniform 
behaviour of the strain distribution near the interface. 41 
At the end we summarize the effects of SI on the en- 
ergy gap by reporting, in Table II, the HOMO-LUMO 
gap-shift resulting from the comparison between the 
suspended-relaxed and -strained NCs, with the corre- 
sponding 17 value. As expected, f2 and the gap-shift 
present the same trend. 





Si 35 -OH 


Sii ( ,-OH 


Si 32 -OH 


Sii 7 -OH 


gap-shift (eV) 


0.07 


0.37 


0.48 


1.29 


Q 


1.2 


1.6 


2.8 


3.0 



TABLE II: Oxydation degree, fi, and HOMO-LUMO gap- 
shift resulting from the comparison between the suspended- 
relaxed and -strained NCs. 



IV. OXIDATION AND STRAIN 

In this section we explore the possibility that the 
strong discontinuity in the energy gap values for the 
NCs analyzed above, could be related to the different 
passivation regime of the interface silicon atoms. It is 
worth noting that a recent size-dependent experimental 
study of Si 2p core-level shift at Si-NC/Si02 interface 



showed that the shell region around the Si-NC bordered 
by Si02 consists of the three Si suboxide states, Si 1+ , 
Si 2+ , and Si 3+ , where for NCs smaller than about 3 
nm the presence of suboxides with higher oxidation 
becomes favoured. 42 In order to focus on this aspect, we 
produced a set of suspended NCs, each with different 
geometry, size, and oxidation degree (see Fig. 4). 
The strain has been intentionally removed from all 
the NCs, through geometry optimization, to elucidate 
further the correlation between the final gap value and 
ft. In parallel to the hydroxided NCs we also built 
the corresponding hydrogenated systems, and we then 
calculated the HOMO-LUMO gaps for all the structures 
(see Fig. 5). What emerges, for the hydroxided case, is 
a strong and clear correlation between Q and the gap 
value. For these sizes, the dependence of the gap on fl 
is much stronger than that due to QC, so that a small 
variation in size could correspond to a large variation 
of the gap value. The effects of QC arc still alive, but 
strongly modulated by that of oxidation. As already 
anticipated, the hydrogenated case does not instead 
depend on il but simply follows the QC. As suggested 
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FIG. 5: (color online) HOMO-LUMO gaps for the hydroxided 
(circles) and hydrogenated (triangles) NCs, together with the 
oxidation/hydrogenation degree Q (squares). 



by other works 23,38 , we expect that for large clusters 
the gap becomes independent on the passivation regime, 
just following the QC. 

Considering the results from the previous Section, we 
could guess that the inclusion of the strain would reduce 
the HOMO-LUMO gaps of Fig. 5, depending (non 
linearly) on the fi value. In our opinion, this would 
reduce the gap-fluctuations amplitude, and also influence 
the QC trend. In fact, by looking at Table I we see 
that, by comparing for example the Sii7 and the Si32 
embedded NCs, the final gap do not follow the QC rule, 
even if they present similar oxydation degrees. In this 
case the strain completely overrules the oxidation effect. 
Instead, for very different Q values the final gap still 
present a clear correlation to the oxidation degree, again 
undependently on the NC size. These considerations 
suggest a final picture in which the gap of small embed- 
ded NCs show an oscillating behaviour with size, and not 
strictly following the QC, as already observed in other 
works. 34,38 The large broadening of the PL spectra, 
observed also at low temperatures, 3,5, 6,8,9,12,43-48 even 
for apparently monodispcrsed multilayered samples, 49 
could be a possible consequence of such fluctuations of 
the electronic gap with the NC size. 

The case of interface defects has been deeply investi- 
gated, sometimes suggesting that they have a small in- 
fluence on the radiative recombination, 9,50-52 and other 
times being indicated as the primary source for radia- 
tive emission. 4,7,10,12,46,47 Moreover, the presence of other 
oxygen-bond types is frequently reported in experimen- 
tal samples, beside the single-bond type discussed here. 
Several works report the Si-O-Si bridge bond as the most 
stable configuration, and the Si=0 double bond as the 
less energetically favoured. 53-56 Besides, these alterna- 
tive bond types are especially favoured for large clus- 
ters sizes. 57,58 In conclusion, we believe that these addi- 
tional parameters (defects, bond types) would not modify 
the general picture related to the parameter Q discussed 
above, making it a crucial quantity for the characteriza- 
tion of small Si-NCs. 



which arc embedded in a silica matrix. The relaxed 
structures show that the surrounding matrix always 
produces some strain on the nanocrystal, especially 
at the Si/Si02 interface. By including the strain on 
the suspended hydroxided nanocrystals we were able 
to reproduce the electronic and optical properties of 
the full Si/Si02 systems. Besides, the removal of such 
strain by a total relaxation will produce a significant 
reduction of the electronic gap. We found evidences 
that the amount of strain exerted on the nanocrystal is 
connected to the interface configuration, in particular to 
the number of oxygens per interface silicons (oxidation 
degree). Moreover we revealed that, while the quantum 
confinement dominates in the hydrogenated nanocrystals 
of all sizes, the behaviour of hydroxided and embedded 
nanocrystals show a strong correlation between the final 
value of the electronic gap and the oxidation degree, that 
seems to override the effects of quantum confinement 
for nanocrystals of diameters below 2 nm. The final 
picture indicates a competition between the oxidation 
and the strain, in which the former tends to blue-shift 
the absorption spectra, and the latter to red-shift it. 
At the end we suggest a final picture in which the gap 
of small embedded NCs shows an oscillating behaviour 
with size, not strictly following the QC rule. The strong 
fluctuations of the electronic gap with the nanocrystal 
size could be responsible for the large broadening of the 
PL spectra that have been observed experimentally, also 
at low temperatures. 3,5,6,8,9,12,43-49 



ACKNOWLEDGEMENTS 



V. CONCLUSIONS 

Considering the /3-cristobalite phase of the Si02 we 
have generated silicon nanocrystals of different size 



This work is supported by PRIN2007 and MAE-CNR 
Italia- Turchia. We acknowledge CINECA CPU time 
granted by INFM (Progetto Calcolo Parallclo). 



1 O. Bisi, S. Ossicini, L. Pavesi, Surf. Sci. Rep. 38, 5 (2000). 

2 L. Pavesi and D.J. Lockwood, Eds., Silicon Photonics, 
Springer, Berlin, 2004. 

3 M. L. Brongersma, P. G. Kik, A. Polman, Appl. Phys. 
Lett. 76, 351-353 (2000). 

4 T. Shimitsu-Iwayama, T. Hama, D. E. Hole, I. W. Boyd, 
Solid-State Electronics 45, 1487-1494 (2001). 

5 A. V. Hernandez, T. V. Torchynska, Y. Matsumoto, S. 



J. Sandoval, M. Dybiec, S. Ostapenko, L. V. Shcherbina, 
Microelectronics Journal 36, 510-513 (2005). 

6 R. Rolver, M. Forst, O. Winkler, B. Spangenberg, H. Kurz, 
J. Vac. Sci. Technol. A 24, 141-145 (2006). 

7 S. Godefroo, M. Hayne, M. Jivanescu, A. Stesmans, 
M. Zacharias, O. I. Lebedev, G. V. Tendeloo, V. V. 
Moshchalkov, Nature Nanotech. 3, 174 (2008). 

8 S. Gardelis, A. G. Nassiopoulou, N. Vouroutzis, N. Frangis, 



5 



J. Appl. Phys. 105, 113509 (2009). 
9 C. Delerue, G. Allan, C. Reynaud, O. Guillois, G. Ledoux, 
F. Huisken, Phys. Rev. B 73, 235318 (2006). 

10 Y. Kanemitsu, Thin Solid Films 276, 44-46 (1996). 

11 I. Sychugov, R. Juhasz, J. Valenta, J. Linnros, Phys. Rev. 
Lett. 94, 087405 (2005). 

12 B. Averboukh, R. Huber, K. W. Cheah, Y. R. Shen, G. G. 
Qin, Z. C. Ma, W. H. Zong, J. Appl. Phys 92, 3564-3568 
(2002). 

13 I. V. Antonova, M. Gulyaev, E. Savir, J. Jedrzejewsky, I. 
Balberg, Phys. Rev. B 77 125318 (2009). 

14 N. Salivati, N. Shuall, E. Baskin, V. Garber, J. M. Mc- 
Crate, J. G. Ekerdt, J. Appl. Phys. 106, 063121 (2009). 

15 R. W. Liptak, U. Kortshagen, S. A. Campbell, J. Appl. 
Phys. 106, 064313 (2009). 

16 P. R. Coxon, Y. Chao, B. R. Horrocks, M. Gass, U. 
Bangert, L. Siller, J. Appl. Phys. 104, 084318 (2008). 

17 D. Konig, J. Rudd, M. A. Green, G. Conibeer, Phys. Rev. 
B 78, 035339 (2008). 

18 O. Lehtonen, D. Sundholm, Phys. Rev. B 74, 045433 
(2006). 

19 L. Sirleto, M. A. Ferrara, I. Rendina, S. N. Basu, J. Warga, 
R. Li, L. Dal Negro, Appl. Phys. Lett. 93, 251104 (2008). 

20 S. A. Empedocles, D. J. Norris, M. G. Bawendi, Phys. Rev. 
Lett. 77, 3873 (1996). 

21 G. M. Credo, M. D. Mason, S. K. Buratto, Appl. Phys. 
Lett. 74, 1978 (1999). 

22 R. Guerra, I. Marri, R. Magri, L. Martin-Samos, O. Pulci, 

E. Degoli, S. Ossicini, Phys. Rev. B 79, 155320 (2009). 

23 I. Vasiliev, J. R. Chelikowsky, R. M. Martin, Phys. Rev. B 
65, 121302(R) (2002). 

24 G. Allan, C. Delerue, M. Lannoo, Phys. Rev. Lett. 76, 
2961 (1996). 

25 H. Kageshima, K. Shiraishi, in: M. Scheffler, R. Zimmer- 
mann (Eds.), Proc. 23rd Int. Conf. Phys. Semicon., World 
Scientific, Singapore, p. 903, (1996). 

26 M. Luppi, S. Ossicini, J. Appl. Phys., 94, 2130 (2003). 

27 M. Luppi, S. Ossicini, Phys. Rev. B 71, 035340 (2005). 

28 F. Djurabekova, K. Nordlund, Phys. Rev. B 77, 115325 
(2008). 

29 P. Ordejon, E. Artacho, J. M. Soler, Phys. Rev. B (Rapid 
Comm.) 53, R10441 (1996). 

30 J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Junquera, 
P. Ordejon, D. Sanchez-Portal, J. Phys.: Condens. Matt. 
14, 2745 (2002). 

31 P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, 
C. Cavazzoni, D. Ceresoli, G. L. Chiarotti, M. Cococ- 
cioni, I. Dabo, A. Dal Corso, S. Fabris, G. Fratesi, S. 
de Gironcoli, R. Gebauer, U. Gerstmann, C. Gougous- 
sis, A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, 

F. Mauri, R. Mazzarello, S. Paolini, A. Pasquarello, L. 
Paulatto, C. Sbraccia, S. Scandolo, G. Sclauzero, A. P. 
Seitsonen, A. Smogunov, P. Umari, R. M. Wentzcov- 



itch, J. Phys.: Condens. Matt. 21, 395502 (2009) http: 
//dx . doi . org/10 . 1088/0953-8984/21/39/395502. 

32 D. E. Yilmaz, C. Bulutay, T. Qagin, Phys. Rev. B 77, 
155306 (2008). 

33 P. Kroll, H. J. Schulte, Phys. Stat. Sol. B 243 (2006). 

34 L. Koponen, L. O. Tunturivuori, M. J. Puska, R. M. Niem- 
inen, Phys. Rev. B 79, 235332 (2009). 

35 K. Seino, F. Bechstedt, P. Kroll, Nanotechnology 20, 
135702 (2009). 

36 X.-H. Peng, S. Ganti, A. Alizadeh, P. Sharma, S. K. Ku- 
mar, S. K. Nayak, Phys. Rev. B 74, 035339 (2006). 

37 A. Franceschetti, Phys. Rev. B 76, 161301 (R) (2007). 

38 Z. Zhou, L. Brus, R. Friesner, Nano Lett. 3, 163 (2003). 

39 M. V. Wolkin, J. Jorne, P. M. Fauchet, G. Allan, C. 
Delerue, Phys. Rev. Lett. 82, 000197 (1999). 

40 K. Dohnalova, K. Kusova, I. Pelant, Appl. Phys. Lett. 94, 
211903 (2009). 

41 D. E. Yilmaz, C. Bulutay, T. Qagin, Appl. Phys. Lett. 94, 
191914 (2009). 

42 S. Kim, M. C. Kim, S-H. Choi, K. J. Kim, H. N. Hwang, 

C. C. Hwang, Appl. Phys. Lett. 91, 103113 (2007). 

43 M. Dovrat, Y. Shalibo, N. Arad, I. Popov, S.-T. Lee, A. 
Sa'ar, Phys. Rev. B 79, 125306 (2009). 

44 M. Sykora, L. Mangolini, R. D. Schaller, U. Kortshagen, 

D. Jurbergs, V. I. Klimov, Phys. Rev. Lett. 100, 067401 
(2008). 

45 S.-W. Lin, D.-H. Chen, Small 5, 72-76 (2009). 

46 Y. Kanemitsu, Y. Fukunishi, M. Iiboshi, S. Okamoto, T. 
Kushida, Physica E 7, 456-460 (2000). 

47 M. Ray, K. Jana, N. R. Bandyopadhyay, S. M. Hossain, D. 
Navarro-Urrios, P. P. Chattyopadhyay, M. A. Green, Solid 
State Communications 149, 352-356 (2009). 

48 G. Belomoin, J. Therrien, M. Nayfeh, Appl. Phys. Lett. 
77, 779-781 (2000). 

49 M. Cazzanelli, D. Navarro-Urrios, F. Riboli, et al., J. Appl. 
Phys. 96, 3164-3171 (2004). 

50 L. E. Ramos, J. Furthmiiller, F. Bechstedt, Phys. Rev. B 
71, 035328 (2005). 

51 J. Derr, K. Dunn, D. Riabinina, F. Martin, M. Chacker, 
F. Rosei, Physica E 41 668-670 (2009). 

52 N. A. Hill, K. B. Whaley, Phys. Rev. Lett. 75, 1130 (1995). 

53 N. B. Nguyen, C. Dufour, S. Petit, J. Phys. Cond. Matt. 
20, 455209 (2008). 

54 D. Jurbergs, E. Rogojina, L. Mangolini, U. Kortshagen, 
Appl. Phys. Lett. 88, 233116 (2006). 

55 K. Seino, J.-M. Wagner, F. Bechstedt, Appl. Phys. Lett. 
90, 253109 (2007). 

56 A. Puzder, A. J. Williamson, J. C. Grossman, G. Galli, J. 
Am. Chem. Soc. 125, 2786-2791 (2003). 

57 G. Hadjisavvas, PC. Kelires, Physica E 38, 99 (2007). 

58 M. Ippolito, S. Meloni, L. Colombo, Appl. Phys. Lett. 93, 
153109 (2008). 



G 



